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ABSTRACT 

Gas-rich primordial disks and tenuous gas-poor debris disks are usually considered as two distinct 
evolutionary phases of the circumstellar matter. Interestingly, the debris disk around the young main- 
sequence star 49 Ceti possesses a substantial amount of molecular gas, and possibly represents the 
missing link between the two phases. Motivated to understand the evolution of the gas component 
in circumstellar disks via finding more 49 Ceti-likc systems, we carried out a CO J— 3— 2 survey 
with Atacama Pathfinder Experiment, targeting 20 infrared-luminous debris disks. These systems fill 
the gap between primordial and old tenuous debris disks in terms of fractional luminosity. Here we 
report on the discovery of a second 49 Ceti-like disk around the 30 Myr old A3-type star HD21997, a 
member of the Columba Association. This system was also detected in the C0(2— 1) transition, and 
the reliable age determination makes it an even clearer example of an old gas-bearing disk than 49 
Ceti. While the fractional luminosities of IID21997 and 49 Ceti are not particularly high, these objects 
seem to harbor the most extended disks within our sample. The double-peaked profiles of HD21997 
were reproduced by a Keplerian disk model combined with the LIME radiative transfer code. Based 
on their similarities, 49 Ceti and IID21997 may be the first representatives of a so far undefined new 
class of relatively old (>8Myr), gaseous dust disks. From our results, neither primordial origin nor 
steady secondary production from icy planetesimals can unequivocally explain the presence of CO gas 
in the disk of HD21997. 

Subject headings: circumstellar matter — infrared: stars — stars: individual (HD21997) 



1. INTRODUCTION 

Most young stars are surrounded by circumstellar 
disks, the natural by-product of star formation. After 
a protostar has formed, its disk plays a crucial role in 
the evolution of the system, first by serving as reser- 
voir for mass accretion, and later by becoming the birth- 
place of the planetary system. At early evolutionary 
stage the mass of the primordial disk is dominated by 
gas, with a few percent of mass in small dust grains. 
The gaseous component plays an import ant role in con- 
trolli ng the dust dynamics in the disk (jBeckwith et al.l 
|2000| ). and its dispersal determines the time available to 
form Jupiter- and Saturn-like giant planets. As the disk 
evolves, the gas content is rem oved by viscous accretion 
(|Lvnden-Bell fc Pringlelll974D a nd /or by photoevapora- 
tion (e.g.. [Alexander et a" 

Current observational results imply that the pri- 
mordial gas becomes la r gely d epleted in th e first 
1 Myr dZ uckcrma n et~all 119951 : iPascucci et"all 120061: 
iFedele et al..,2010, ). A gas-poor disk appears, where the 
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lifetime of individual dust grains under the influence of 
radiative forces - without the stabilization effect of the 
surrounding gas - is far less than the stellar age. The 
dust grains are believed to be continuously replenished by 
collisions and/or evapo ration of previously formed plan- 
etesimals (|Wvattll2008l and references therein). In these 
debris disks only a small amount of gas is expected. Sim- 
ilarly to the dust, this gas could be second generation, 
produced by sublimati on of p lanetesimals, photodcsorp- 
tion from dust grains (iGrigor ieva e t al. 2007) , or vapor- 
ization of colliding dust particles ( Czechowski fc Mamil 
|2003). 

So far only a few debris disks are known with de- 
tectable gas component. The edge-on orientation of the 
disk around /?Pic allowed the detection of a very small 
amount of circumstellar gas (iVco ~ 6 x lO^^cm"^) 
through the presence of absorption/emi ssi on lines 
•Roberge et all l2000l: [BTandeker et afl [200l. IRedfieldl 



20071 ) successfully exploited the favorable edge-on ge- 



ometry of the disk around IID32297 to detect gas via 
Na I. In contrast to the disks mentioned above, the de- 
bris disk around the young main-sequence star 49 Ceti 
seems to have substantial (^13 Mm) molecular gas 
(IZuckerman et all 119951 : IDent et all 120051: IHughes et all 
I2008D. The" origin of the gas in the above mentioned 
systems is currently under debate. It can be residual 
primordial gas that s urvived longer in the outer disks 
than usually assumed (jKrivov et al.1 12009') or it may have 
formed or bee n released recent l y. Ba sed on dynami- 
cal arguments, [Fernandez et al.l (|2006[ ) suggested that 
the gas in f3 Pic is secondary. Analyzing high-resolution 
data obtained with the SMA interferometer at 230 GHz, 
IHughes et al.l (|2008[ ) proposed that 49 Ceti is in a late 
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stage of its evolution, possibly representing the link be- 
tween gas-rich pri mordial disk s and gas-poor optically 
thin debris disks. iRhee et aL (120071) suggeste d an age 
of 20 Myr for 49 Ceti. while Ixhi et all (|2001h derived 
^8 Myr as the age of the star. In the case of the older age, 
49 Ceti would challenge the current picture of gas disk 
evolution, while the lower value could still be marginally 
consistent with a primordial disk phase. The confirma- 
tion of the existence of debris disks containing a signifi- 
cant amount of gas would require to find and study more 
49 Ceti-like systems with reliable age estimates. It might 
well be possible that a number of similar systems ex- 
ists, since most debris disks that are similar to 49 Ceti in 
terms of age and fractional luminosity have never been 
observed in molecular lines (most such candidates are in 
the southern hemisphere). Motivated by this fact, we 
carried out a survey with the APE>s0 radio telescope to 
detect molecular gas in 20 infrared-luminous debris disks. 
Here we review the results of this survey and report on 
the discovery of a second 49 Ceti-like disk around the 
30Myr-old star HD21997. 

2. SAMPLE DESCRIPTION 

For candidate selection we adopted the 49 Ceti system 
as a template. This Al-type star harbors a dusty disk 
that re-emits a fraction /dust ~ 8 — 9 x lO^*^ of the star's 
radiation at infrared wavelengths. This fractional lumi- 
nosity is an order of magnitude lower than the corre- 
sponding value in primordial disks. We used the follow- 
ing target selection criteria: (1) spectral type between 
AO and F6; (2) /dust in the range of 5 x IQ-Hob x lO'^ 
that excludes both primordial disks and very tenuous 
debris disks; (3) the excess emission is confirmed by an 
instrument independent of IRAS (4) ages between 10 and 
100 Myr, the age estimate is based on kinematic group 
membership or other reliable dating metho ds. In total, 
20 candidates were selected from the lists of iChen et al.l 
11005) and Moor ct al. (2006, 2011). Their basic proper- 
ties are given in Table [ij 

For most sources, disk parameters (temperature, ra- 
dius, and fractional luminosity) were adopted from the 
literature (Table [1]). In all cases, disk radii were esti- 
mated by adopting blackbody-like dust emission. For 
IID121617 no previous literature data were found, thus 
we collected infrared photometry from th e IRAS FSC, 
AKARI IRC, AKARI FIS and WISE (jWright et alj 
[2OT0h . For HD21997, we also compiled a new spectral 
energy distribution based on literature data and our own 
photometry derived from archival observations obtained 
with the Multiband Imaging Photometer for Spitzer 
(see Figure [H 55.3±2.2mJy at 24 /im and 662±47mJy 
at 70 ^m; for t he description of the processing see 
iMoor et al.ll201ll ) . Disk parameters for the se two targets 
were derived following I Moor et al.l (|2011[ ) . For sources 
where submillimeter observations were available, we com- 
puted dust masses assuming opticall y thin emission with 
i^svo^m = 2 cm^ g~^ and 13 — 1 (Nilsson et al.ll20ldl ). and 
dust temperature from Table [1] For comparison of the 
fundamental parameters, 49 Ceti is also added to Table[TJ 

® This publication is based on data acquired with the Ata- 
cama Pathfinder Experiment (APEX). APEX is a collaboration 
between the Max-Planck-Institut fiir Radioastronomie, the Euro- 
pean Southern Observatory, and the Onsala Space Observatory. 



3. OBSERVATIONS AND DATA REDUCTION 

Our su rvey was carr ied out with the APEX 12 m 
telescope (jGiisten et al.l '20061 in service mode, between 
2008 October and 2009 November. All objects were 
observed at the 345.796GH z ^'^CO J=3-2 line using 
the SHeFI/APEX2 receiver (jVassilev et al.l [20081 ). One 
source, HD21997 was also observed in the J=2— 1 tran- 
sition of i^CO (at 230.538 GHz) with the SHeFI/APEXl 
receiver. For the backend, we used the Fast Fourier 
Transform Spectrometer with 2048 channels providing a 
velocity resolution of 0.42 and 0.64 kms^^ in the J=3— 2 
and J=2— 1 transitions, respectively. An on-off observ- 
ing pattern was utilized with beam switching. The total 
on-source integration time for most sources ranged be- 
tween 10 and 30 minutes. For HD21997, we integrated 
longer for a better characterization of the line profile. 
The weather conditions were generally dry, the precip- 
itable water vapor was below 1.3 mm for most of the 
J=3— 2 observations and ranged between 1.3 and 2.7mm 
during the J=2— 1 measurements. The data reduction 
was performed using the GILD AS / CLASS package^- For 
the final average spectrum, we omitted noisy scans and a 
linear baseline was subtracted from each individual scan. 

4. RESULTS AND ANALYSIS 

Among the 20 targets, one system, HD21997, was de- 
tected at > 5(7 level in both CO lines. Figure [5] shows the 
baseline-corrected CO profiles for HD21997. Both lines 
display double peaked line profile with identical peak po- 
sitions. The central velocities of both lines are consistent 
with the systemic veloci ty of the star (-f 17.3±0.8kms^^; 
iKharchenko et all l2007| ). We integrated the intensi- 
ties/fluxes over an interval of 8kms~^ that covers the 
whole line profile. The beam efficiencies and Kelvin-to- 
Jansky conversion factors were taken from the APEX 
web pagj*^. For the non-detected sources, upper lim- 
its were estimated as TimsAv-v/iV, where Trms is the 
measured noise, Av is the velocity channel width, and 
TV is the number of velocity channels over an interval 
of lOkms"^. The total mass (or upper limit) of CO 
molecules (Afco) was estimated assuming optically thin 
emission and local thermodynamic equilibrium (LTE). 
The excitation temperature (Tex) was assumed to be 
equal to the dust temperature in Table [1] (i.e., gas and 
dust are sufficiently coupled). The obtained line intensi- 
ties/fluxes as well as the estimated CO masses are listed 
in Tabled Note that for HD21997 the CO masses com- 
puted independently from the (2—1) and (3—2) transi- 
tions are significantly different. 

Figure [3Uleft panel) shows the integrated C0(3— 2) 
fiuxes and upper limits, normalized to 100 pc, plot- 
ted against the fractional luminosities of the disks. 
For comparison, additional protoplanetary /debris disks 
around A0-F6-type pre- main /main- s equence sta r s, in- 
cluding 49 Ceti, are also displayed (jDent et al.l 120051 : 
IGreaves et al.ll20d0airbl ) . Our observations fill the gap be- 
tween Herbig Ae/Be disks and older debris disks. Note 
that the fractional luminosities of HD21997 and 49 Ceti 
are modest even within the debris sample, and signif- 
icantly lower than those of the primordial disks. Thus, 
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/dust does not appear to be a good proxy for the presence 
of CO gas in debris disks. Figure [Sjright panel) presents 
the integrated C0(3— 2) fluxes versus the (blackbody) 
radii of the dust disks. Interestingly, the two definite de- 
tections, HD21997 and 49 Ceti, harbor the most extended 
disks, suggesting that large radius and low dust temper- 
ature may be essential for CO detection. Although the 
radii in this analysis rely on the assumption of blackbody 
grains, the conclusion holds in the case of realistic grain 
size distributions. Using the blackbody assumption, we 
may systematically underestimate the true radius, due 
to the presence of inefficiently emitting small particles. 
However, assuming similar grains in all disks, the rela- 
tive distribution of disk radii would not differ from the 
blackbody case (see e.g.. lWvattll2008D . 

HD2 1997 is an A3-ty pe star at a distance of 
72 pc (jvan LeeuwenI l2007fl. a meniber of the ^30Myr 
old Columba group (|Mo6r et"ari l2006t iTorres et al.l 
I2008D. Fitting an ATLAS9 atmosphere model 
(|Castelli fc Kuruca I2003D to the optical and near-IR 
{Hipparcos, Tycho-2, Two Micron All Sky Survey) data, 
assuming solar metallicity and log g = 4.25 yields T pff = 
8300 K. The evolutionary tracks of S iess et all (|2000f ) im- 
ply a stellar mass of 1.85 Mq. We modeled the measured 
line profiles of HD21997 with a simple disk geometry as- 
suming a combination of a radial power-law and a vertical 
Gaussian profile for the density distribution: 

ncoir, z) = nco,in(r/i?i„)-"e-^'/2^' . (1) 

We fixed the following parameters: H — Q.\r, a = —2.5, 
i?i„ = 63 AU (Table ID, and Rr,ut = 200 AU (ty pical for 
Herbig Ae/Be disks: (Panic fc Hogerheiid3l2009[ ). We as- 
sumed an H2 abundance relative to CO of 10"*^, and that 
gas and dust grains - the latter act like blackbody - are 
well mixed, prescribing that the gas kinetic temperature 
and dust temperature distributions are identical. The 
velocity of the material in the disk was derived by as- 
suming Keplerian orbits around a star of 1.85 Mq mass. 
Then the CO level populations at each disk position, and 
the resulting emission line profiles were calculated using 
the non-LTE spectral line ra diation transfer code LIME 
(|Brinch fc Hogerheiidell2010D . First, we fitted the (3-2) 
line by adjusting nco,in and the disk inclination. The 
best-fitting model spectrum with nco.in = 10 cm~'^ and 
i = 45° is overplotted with dashed line in Figure [2j With 
the same parameters we also computed a CO (2-1) pro- 
file. As Figure [5] shows, this model significantly under- 
estimates the observed C0(2-l) feature. The reason for 
this discrepancy is the same as what causes the difference 
in the CO mass estimates from (2-1) and (3-2) lines: 
the ratio of integrated C0(3— 2) flux to the integrated 
C0(2— 1) flux is only 1.43±0.37, significantly lower than 
the ratio of 3.8, expected for Tex ~ 60K in LTE condi- 
tion. This low line fiux ratio corresponds to an excita- 
tion temperature of 13.1±2.7K. This is unrealistically 
low (subthermal) for being an LTE value, suggesting 
that the density of collision partners (H2) is lower than 
the critical density, and the excitation is not collisionally 
dominated. In order to provide a model that can fit both 
lines simultaneously, we gradually decreased the H2/CO 
abundance ratio and repeated the above-mentioned mod- 
eling process. We found that with H2/CO=1000±500, 
'T'CO.in = 22 ± 5cm~'^, and i — 45lJoi both line profiles 



can be fitted (solid line in Figured]). Note that in this 
non-LTE model the kinetic and excitation temperatures 
of the gas are different. The total CO mass predicted by 
this model is Mqo = 3.5 x lO-'^M®. 

5. DISCUSSION 

Its reliable age determination makes IID21997 the old- 
est example for a gas-bearing debris disk. In many as- 
pects it resembles the somewhat younger 49 Ceti. Both 
systems contain an A-type central star that produces en- 
ergetic UV that could dissociate CO molecules in the 
vicinity of the star. 49 Ceti and HD21997 clearly stand 
out from our sample in terms of disk radius, and of har- 
boring a large amount of relatively cold dust (Tdust < 
80 K and Mdust ~ 0.1 M^). Note also that these two 
systems exhibit very similar Mco/-^dust ratios (^0.003). 
Based on their similarities, we speculate that HD21997 
and 49 Ceti may be the first representatives of a so far 
undefined new class of relatively old (>8Myr), gaseous 
d ust disks. 

iHughes et al.l ()2008[ ) claim that the disk around 49 Ceti 
contains predominantly primordial gas and may repre- 
sent the link between gas-rich primordial disks and gas- 
poor debris disks. It is a question whether HD21997 
system may be of similar origin. The gas clearing pro- 
cess in primordial disks is expected to progress outwards 
due to photoevaporatio n driven by th e central star (e.g., 
Alexander et al. 2006; P ascucci fc Tac hibana 2010), thus 
the last reservoir of gas will be the outermost part of the 
circumstellar disk. Indeed, IID21997 and 49 Ceti possess 
the largest disks and consequently have the longest ex- 
pected survival time for gas in our sample (though the 
evaporation timescale also depends on the high energy 
fiux of the central star). A confirmed primordial origin 
of the gas in the IID21997 system, would pose a serious 
question to the current paradigm since its age of ~30 Myr 
significantly exceeds both the model predictions for disk 
clearing and the ages of the oldes t T Tauri-like or tran si- 
tional gas disks in the literature ([Kastner et al.|[2"008f) . 

Primordial CO gas can survive only in the case of ef- 
ficient shielding from the stellar/interstellar high-energy 
photons. We determined the stellar UV fiux from the fit- 
ted ATLAS9 atmosphere model (see Section and the 
UV compone nt of t he interstellar radiation field (ISRF) 
from HabinJ (|1968f ). For each disk position, where the 
II2 and CO column densities are provided by our model, 
we analyzed the shieldi ng effici e ncy u sing the photodis- 
sociation model of Viss er et al.l (|2009( ) for different A'^co 
and A^Ha pairs (shielding by dust grains is negligible 
in this tenuous disk). We found that no region in the 
disk is shielded enough to provide a CO lifetime longer 
than 500 yr. Adopting a lower scale height would lead 
to higher radial column densities but would not affect 
significantly the vertical column densities, thus the UV 
photons of the ISRF - which dominate the stellar UV 
flux almost everywhere in the disk - could efficiently pho- 
todissociate CO molecules. In the course of modeling we 
assumed that the gas and dust are sufficiently coupled 
leading to a common temp erature, but in tenu ous debris 
disks this may not be tr ue (jKamp fc van Zadel hoff 2005 
iZagorovskv et al1l2010D . Assuming a lower gas temper- 
ature similar to the measured excitation temperature 
would allow the existence of a larger amount of hydro- 
gen gas. However, the gas is not likely to cool down to 
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such a low temperature all over the disk. Thus based 
on this result, as well as on the obtained H2/CO ratio 
of '--^lOOO that is lower than expected for a primordial 
composition, the scenario of primordial origin of gas in 
HD21997 is unlikely. 

Is it possible then, that the gas in this disk is of sec- 
ondary origin, being continuously replenished from icy 
planetesimals? In this scenario the gas may have a 
very low H2/CO ratio. Without the presence of a large 
amount of H2, shielding against UV photons is weak. 
Our modeling predicts CO photodissociation timescales 
of less than 500 years in the whole disk. In order to repro- 
duce the observed CO mass of ~3.5x lO"** Mq, CO has to 
be released from solids with a rate of >7x 10"'' M^yr"^. 
Pure CO ice evaporates at temperatures above 18 K, thus 
it is volatile even far from the luminous central star, 
meaning that the surface of planetesimals is very likely 
already depleted in CO ice. However, CO ice can per- 
sist in deeper layers and/or can be trapped in mixed 
H2O-CO ices even on the surface at temperatures be- 
low ^140 K. Destructive collisions between planetesimals 
can lead to the release of subsurface ices. In addition, 
frequent collisions of icy grains with smaller particles - 
mainly with /3 meteoroids - could produce a continu- 
ous flux of CO via vaporization. Photodesorption from 
solids can also significantly contribute to the gas produc- 
tion. Extrapolation of the current production rate for the 
last 20 Myr (assuming a primordial gas-rich disk phase in 
the first 10 Myr and a steady-state disk evolution after- 
wards) would yield a total of >14M0 of CO released 
from planetesimals/grains. Adopting a CO mass abun- 
dance of 10% in the planetesimals (see the composition 



of the comet IIale-Bopp: lHuebner fc Benkhoflri999D . this 
scenario would require the complete destruction of more 
than 140 mass of planetesimals in the outer disk. It 
would significantly exceed the full initial dust content of 
a typical protoplanetary disk, making this steady-state 
scenario questionable. A more satisfactory explanation 
would be that the system is currently undergoing a phase 
of temporarily high CO production. The origin of this 
contemporeous gas production might be imprinted in the 
spatial distribution of the gas and dust, and could be re- 
vealed with future interferometers. 

Our results indicate that neither primordial origin nor 
steady secondary production can unequivocally explain 
the presence of CO gas in the disk of HD21997. An 
on-going temporarily high CO production may be more 
likely. Detection of other gas components and transi- 
tions with Herschel and ALMA, as well as the better 
characterization of the disk structure may lead to the 
deeper understanding of this enigmatic system and clar- 
ify whether 49Ceti and HD21997 are the first examples 
of a so far less studied phase of disk evolution. 
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TABLE 1 
Stellar/disk properties. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(S) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


ID 


Scl. 


SpT 


D 


Mcmb. 


^rad 


^dust 


-f^dust 


/dust 


Rcf. 


^^du.t 


Rcf. 


ton 


ico , 


Sco 










(pc) 




(kms-1) 


(K) 


(AU) 


(10-4) 




(M©) 




(minute) 


(Kkms-l) 


(Jykms-1) 


(lO-^M^) 


HD 3670 


1 


F5V 


(76) 


ColA 


+ S.6 


53 


41 


5.4 


1 






35.5 


<0.050 


<1.50 


<0.85 


HD 15115 


1 


F2 


45 


BPMG 


+ 8.8 


56 


42 


4.8 


1 


0.036 


1 


39.4 


<0.045 


<1.35 


<0.28 


HD 21997 


2 


A3IV/V 


72 


ColA 


+ 17.3 


64 


63 


5.9 


5 


0.14 


1 


255.3 


0.109 + 0.014 


3.28 + 0.53 


1.79 + 0.32 


HD 21997'* 
























115.9 


0.078 + 0.014 


2.29 + 0.47 


4.9+1.1 


HD 30447*^ 


1 


F3V 


SO 


ColA 


+ 21.3 


62 


38 


8.8 


1 


0.16 


2 


24.8 


<0.050 


<1.50 


<1.01 


HD 32297^ 


2 


AOV 


112 




+ 23.0 


85 


28 


54.0 


2 


0.55 


3 


73.6 


<0.027 


<0.80 


<1.25 


HD 35841 


1 


F5V 


(96) 


ColA 


+ 23.1 


68 


23 


15.2 


1 






13.5 


<0.065 


<1.95 


<1.96 


HD 38207 






(93) 














<0 28 












HD 106906 


3 


F5V 


92 


LCC 


+ 11.1 


90 


20 


14.0 


4 






24.2 


<0.069 


<2.08 


<2.24 


HD 110058 


2 


AOV 


107 


LCC 


+ 5.0 


130 


11 


25.0 


2 


<0.34 


5 


43.0 


<0.051 


<1.54 


<2.90 


HD 113766 


3 


F4V 


123 


LCC 


-0.6 


330 


3 


21.0 


4 






28.0 


<0.055 


<1.64 


<8.75 


HD 114082 


3 


F3V 


85 


LCC 


+ 5.2 


110 


10 


30.0 


4 






20.6 


<0.066 


<1.98 


<2.11 


HD 115600 


3 


F2IV/V 


110 


LCC 


+ 4.7 


120 


10 


16.0 


4 






24.0 


<0.056 


<1.67 


<3.15 


HD 117214 


3 


F6V 


110 


LCC 


+ 7.2 


110 


11*= 


7.0 


4 






10.8 


<0.078 


<2.33 


<4.12 


HD 121617 


2 


AlV 


(120) 


UCL 


+ 13.6 


105 


28 


48.0 


5 


<0.33 


6 


14.2 


<0.070 


<2.09 


<4.24 


HD 164249 


2 


F5V 


48 


BPMG 


-0.2 


70 


27 


10.0 


2 


<0.076 


7 


6.7 


<0.129 


<3.85 


<0.99 


HD 172555 


2 


A7V 


29 


BPMG 


+ 2.0 


320 


2 


8.1 


2 






26.2 


<0.047 


<1.39 


<0.39 


HD 181327 


2 


F6V 


52 


BPMG 


+ 0.2 


75 


25 


35.0 


2 


0.40 


7 


31.6 


<0.037 


<1.11 


<0.34 


HD 191089^ 


2 


F5V 


52 


BPMG 


-5.8 


95 


15 


14.0 


2 


0.022 


4 


17.1 


<0.085 


<2.55 


<0.92 


HD 192758 


1 


FOV 


(62) 


Argus 


-11.1 


56 


52 


5.4 


1 






9.8 


<0.093 


<2.78 


<1.07 


HD 221853^ 


1 


FO 


68 


LA 


-4.2 


83 


22 


7.9 


1 






3.4 


<0.131 


<3.91 


<2.24 


49 Ceti 




AlV 


59 






80 


60 


7.9 


2 


0.074 


8 








2.60±0.54'^ 



Note. — Column 1: identification. Column 2: reference for the selection. l 4Moor et a.l.|f^7mi) . 2- FCT^or et al.H^OQgl). 34Clien et al.M2Q0S|) - Column 3: spectral type. Column 4: 
distance. Parenthesis indicate photometric or kinematic distances, otherwise Hipparcos distances from [vanLeeuweT] 20U71 ) areused. Column 5: membership status. ColA: 
Columba Association; BPMG: /3 Pic moving group, LCC: Lower Centaurus Crux association, UCL: Upper Centaurus Lupus association, Argus: Argus moving group, LA: Local 
Association. Co lumn 6: heliocentric radial ve locity . Colu mn 7: disk temperature. Coiu mn8: disk radius . Colu mn9: fractional dust luminosity. ColumnlO; references for disk 
parameters. (l) |Moor 5 t al. |2011), (2) Rhce e t al.| HMTD , (3) millenbrand et ^ ji^OO^j). (4)|Chen et al.| ^jjOOSI). (5) this work . Co lumn 11: dust mass. Col umn 12: reference for 
submi llime ter measurem ent that was used in t he dust mass esti mate. ( 1) 1 WiiliamsU; Andrewj j^iji^^) . (2)p'ilsson et aLinTTTD , (3') |Maness et al.| nHH^ ) , (4) |Hoccatagliata et~^ 
<^005|) . (5) ISylvester et al.j 42001i) . (6) iSherct et al.| r5004|) . (7) INilsson et al.| (8) |jjong e^dl (^U^ ) . Column 13: on-source integration time. Column 14: line intensity 

for CO J=3 — 2. Intensity units are main-beam brightness temperature. Column 15: integrated line flux. Column 16: estimated mass of the CO gas. 
^ Parameters from the J=ii2-1 line. 

^ IT^astner et al.| f^01l!l|) reported the nondetection of CO(2-l) emission. 

For those objects where our distance estimate differs from the literature value by >10%, the Rjust rcscalcd according to our distance. 

^ CO mass of 49 Ceti was derived using the CO(3-2) observation of |Dent et al.| <500S| ). 
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Fig. 1. — Spectral energy distribution of HD21997, Photometric data presented in this Ggure are color-corrected. 
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Fig. 2. — CO spectra of HD21997. The CO(3—2) spectrum is binned by a factor of two. Dotted line marks the radial velocity of the star. 
The dashed line (blue in online version) corresponds to a disk model with H2/CO=10^ , the solid line (red in online version) to H2/CO=lCP . 
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Fig. 3. — Left panel: in tegrated CO (3— 2) Buxcs or upper limits for our sample, and for protoplanetary /debris disks around A0-F6-type 
stars from literatu re data IDent et al£ 2005: Greaves ct al. 2000a b) normalized to lOOpc arc plotted against fractional luminosities. 49 Ceti 
fDent et al\\200'B) . and HD21997 are plotted with larger filled symbols. Right panel: integrated CO(3—2) fluxes for our sample and for 
49 Ceti as a function of disk radii. 



